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ABSTRACT 

Using data from the Sloan Digital Sky Survey data release 3 (SDSS DR3) we investigate 
how narrow (<700km/s) C IV and Mg II quasar absorption line systems are distributed around 
quasars. The C IV absorbers lie in the redshift range 1.6 < z < 4 and the Mg II absorbers in 
the range 0.4 < z < 2.2. By correlating absorbers with quasars on different but neighbouring 
lines-of-sight, we measure the clustering of absorbers around quasars on comoving scales be- 
tween 4 and 30Mpc. The observed comoving correlation lengths are Tq ^ 5/i^^Mpc, similar 
to those observed for bright galaxies at these redshifts. Comparing with correlations between 
absorbers and the quasars in whose spectra they are identified then implies: (i) that quasars 
destroy absorbers to comoving distances of ^ SOOkpc (C IV) and ~ SOOkpc (Mg ll) along 
their lines-of-sight; (ii) that > 40% of C IV absorbers within 3,000 km/s of the QSO are not 
a result of large-scale clustering but rather are directly associated with the quasar itself; (iii) 
that this intrinsic absorber population extends to outflow velocities of order 12,000 km/s; (iv) 
that this outflow component is present in both radio-loud and radio-quiet quasars; and (v) that 
a small high-velocity outflow component is observed in the Mg II population, but any further 
intrinsic absorber component is undetectable in our clustering analysis. We also find an indi- 
cation that absorption systems within 3,000 km/s are more abundant for radio-loud than for 
radio-quiet quasars. This suggests either that radio-loud objects live in more massive halos, or 
that their radio activity generates an additional low-velocity outflow, or both. 

Key words: galaxies:active, quasars:absorption lines; accretion, accretion discs; large-scale 
structure of Universe 



1 INTRODUCTION 

Intervening absorption lines in QSO spectra provide a wealth of in- 
formation on the gaseous Universe from high redshift to the present 
day. They also allow us to probe the metallicity and ionisation state 
of the gas in environments ranging from voids to galaxy halos and 
disks. QSO absorption line systems (QSOALSs) are attributed to 
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two main sources: (i) material associated with the host galaxy of the 
QSO : either in radiatively driven winds close to the accretion disk 
around the black hole or in outflowing material in the host galaxy 
and its surroundings ; (ii) intervening galaxies along the line-of- 
sight. 

QSOALSs are usually split into three classes: (i) systems 
within a few thousand km/s of the QSO systemic redshift, termed 
"associated" systems; (ii) systems with velocity differences in the 
range of 10,000 to 60,000 km/s; (iii) systems beyond this velocity 
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that are generally assumed to be spatially disconnected from the 
QSO. Additionally, QSOALSs are split into two different classes 
depending on their line width. Broad absorption lines (BALs) are 
unambiguously associated with outflowing material from the cen- 
tral region of the QSO; narrow absorption lines (NALs), with typ- 
ical line widths of several hundred km/s, can be associated either 
with the QSO host or with intervening systems. 

It is generally agreed that strong (rest equivalent width 
[EW]>0.3A) metal NALs located at large velocity separations 
from the QSO probe diffu se gas in the halos (and possibly disks) 
of ordinary galaxies (e.g. [ Bergeron & Boiss^ 1 19841 : ISteidel et al.l 
1 19971 ; [Tinker & Ched 120070 . The currently favoured model pos- 
tulates the existence of small (of order 1-lOkpc diameter) clouds 
distributed throughout a galaxy's dark matter halo (see e.g. 
IChelouche etal] l2007h . The best studied transition is that of 
Mg nAA2796.4, 2803.5. The imaging studies of ISteidell j 19931) 
suggested that Mg II extends out to a radius of ~40/i~^ kpc in 
galaxy halos. Recent results show that Mg II absorbing gas extends 
to distances greater than 100 kpc, but with lower covering fr action 



flicting results and conclusions if Weymann etal. 1979; Foltz et alj 
ll986l ; lYoung et alJl982l ; ISargent et al..l988. ; ,Ellingson et al..l994ir 



teowen et ani2006l ; IZibetti et alj2007l ; lKacprzak etai]|2007h . 

The picture becomes much more complicated as we look 
closer towards the QSO in velocity space. A substantial body 
of work has focused on BALs, which cannot be associated with 
intervening galaxies, because they have line widths of at least 
1000 km/s. Commonly reaching sub-relativistic velocities of up 
to 30,000 km/s BALs are u nambiguously associated with AGN- 
driven outflowing gas (e.g. [W evmann et al. 1979; Yuan & Wills! 
l2003l;lRichardsll2006l ; lGangulv et al. 2 007; Lundgren et al, 20071) . 

Many NALs are known to be caused by ordinary galaxies, 
and therefore any additional population associated with the QSO 
must be studied statistically. In general, one looks for an excess 
of NALs relative to the background density in the vicinity of the 
QSO. Because galaxies are clustered, however, one expects an ex- 
cess of NALs due to galaxies surrounding the QSO. This excess 
is expected to extend out to distances of a few galaxy correlation 
lengths (~ 20 Mpc, or velocity separations of ^^3000 km/s). QSO 
absorbers at velocity separations less than this nominal value of 
^ 3000 km/s are termed associated absorption lines (AALs). 

Some authors have dealt with the clustered component by 
restricting their analysis to velocity differences t hat are larger 
than th is value. For example, iRichards et al ]( ll999l) and lRichardj 
( l200lh report an excess of narrow C IV absorbers with veloci- 
ties >5000km/s in radio quiet QSOs as compared to flat-spectrum 
radio-loud QSOs. Other studies have attempted to gain insight 
into the physical nature of AALs in order to ascertain whether 
or not a proportion of these systems may be associated with the 
QSO. Based on line dia gno stics and time variability a rguments, 
iBarlow & SargentI ( 1 19971) and lHamann et alj ( ll99l[200ll) conclude 
that some AALs lie in the QSO host galaxy itself. Lo wer ionisation 
AAL s probably lie at great e r dis tances from the QSO jBarlow et al.l 
1 1997b . IVanden Berk et"al] 1 I2OO&) find that associated Mg II ab- 
sorbers have higher ionisation than intervening systems, clearly 
showing that ionising radiation from the QSO can affect the state 
of the surrounding gas. The fraction of AALs in steep- and flat- 
spectrum radio-loud, and ra dio-quiet QSOs, may also provide clues 
to the origin of the AALs 1 Foltz et al. Il 19881; lAldcroft et al.lll994 
Ganguly et al.l200ll ; lB^er et alj2002l ;lvestergaard 2003; Richards 



200 ih . However, the number of radio-loud QSOs studied thus far 



has been quite small and the results obtained have been inconclu- 
sive. 

There have been many studies of the clustering of NALs 
around small samples of QSOs. These studies have yielded con- 



Recent studies with larger QSO catalogues have generally found 
that t here is an excess NAL population closely associated with the 
QSO ( lRichardsll200ll ; IVestergaar d 2003), but there is still ambigu- 
ity as to whether this population arises from neighbouring galaxies 
or from gas associated with the QSO, its host galaxy and its halo. 

The nature of both AALs and higher velocity associated ab- 
sorbers has potentially important implications for understanding 
feedback processes in QSOs. Feedback is a key ingredient in mod- 
em day galaxy evolution models, invoked to explain observations 
as diverse as the shape of the galaxy luminosity function, the appar- 
ent bimodality in the colour distributions of galaxies, the entropy 
of the intracluster gas, and t he tight correlation between black hol e 
and galaxy bulge mass (e.g. lCroton et alj2006l ; lBower et alj200^ . 
The basic principle is that a large source of energy, originating from 
either an intense burst of star formation or an Active Galactic Nu- 
cleus (AGN), causes the expulsion of gas from the galaxy. This 
feedback halts star formation and heats the surrounding intergalac- 
tic medium (IGM). 

A few studies of both low- and high-redshift radio galaxies 
have yielded direct observational evidence for signifi cant (in terms 
of to t al energy and mass) AGN-driven outflows tEmonts et al.l 
I2OO5I ; iMorganti et ai]|2005l I200I iNesvadba et al.ll2006l.l2007l) , but 
such syst ems are very rar e . BALs are ob served i n 15-20% of 
all (jSOs jFoltz et al.Hl990l ; lHewett & Folt z 2003; R eichard et all 
I2OO3I ; iTrump et al.ll2006l) . but these are best understood in terms 
of winds produced in the immediat e vicinity of the black hole 
jMurrav & Chiang|ll995l ; lElvisll2OO0l) . The impact of such winds 
on the intergalactic medium of the galaxy is not well understood 
(Wevmann 2002). 

Associated narrow absorption lines are tanta lizing, because 
they m ay be the signature of galaxy-wide feedback jTremonti et al.l 
I2OO7I) . However, in order to reveal the true outflow signal, the con- 
tribution from galaxies clustered around the QSO must first be 
subtracted. This can be done if one is able to compute the cross- 
correlation between QSOs and absorbers on neighbouring lines- 
of-sight, since such correlations canno t be due to processes occur- 
ing within the QSO or its host galaxy. lOutram et al.l (|200ll) cross- 
correlated QSOs and metal absorption systems detected in the 2QZ 
QSO survey. They detected a marginally significant positive corre- 
lation on scales of 10/i~^ Mpc. lHennawi & Prochask3 ( l2007h used 
17 Lyman Limit absorption systems with transverse separations 
from the QSO of less than 6 Mpc to measure the clustering of op- 
tically thick neutral hydrogen around QSOs. Their inferred cluster- 
ing amplitude implied a 15-50% underdensity of line-of-sight ab- 
sorption systems, likely caused by photoevaporation of absorbers 
by the ionising flux of the QSO. 

In this paper, we study the clustering properties of a sample of 
6456 C IV systems and 16137 Mg II systems drawn from the SDSS 
Data Release 3 (DR3) QSO catalogue. Our aim is to recover the 
fraction of AALs that are truly intrinsic to the QSOs. We perform 
a 3 dimensional QSO-absorber cross-correlation analysis to mea- 
sure the contribution to the observed NALs from galaxies clustered 
around the QSOs. We compare and contrast our results for Mg II 
and C IV absorber samples. These ions have very different ionisa- 
tion energies, 15.03eV and 64.5eV respectively, and thus provide 
information on the physical state of the gas around QSOs. In Sec- 
tion [3] we describe the absorber selection and the basic properties 
of our samples. Section|4]presents the line-of-sight correlation be- 
tween absorbers and QSOs in redshift space. In Section |5] we es- 
timate the large scale Mg II-QSO and C IV-QSO correlations. In 
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Section |6] these 3-D correlations are used to estimate the excess 
number of absorbers along the QSO line-of-sight due to clustering. 
The main results are presented in Section |7] including the differ- 
ence in NAL distributions for radio-loud and radio-quiet QSOs. We 
conclude and summarise in Sections [H and |9] 

Throughout the paper we assume a flat cosmology with Q,a = 
0.7, f^M = 0.3, Hq = 100/ikms"^ Mpc"^ Unless otherwise 
stated distances are given in comoving units assuming h = 0.7. 



2 THE SLOAN DIGITAL SKY SURVEY 



The S DSS jYork et aljr2000h is using a CCD ca mera dOunn et al.1 
on a 2.5-m telescope dOunn et alll2006l) at Apache Point 
observatory to perform a 5-band jpukugita et alj|l996) photomet- 
ric and spectroscopic survey of the high Galactic latitude sky. 
Photometric calibration is provided by simulta neous observations 
with a 20-inch telescope at the same site (see iHogg et al.|[200ll : 
ISmith et al.|[2002l ; Istoughton et"ai]|2002l : iTucker et al.ll2006h . The 
survey data-processing software measures the properties of each 
detected object in the i maging data, and determines and a pplies 
photometric calibrations jPier et aill2003l : llvezic et alj|2004l) . The 
spectroscopic survey provides medium resolution (R~2000) spec- 
tra with a wavelength cover age f rom 3800 to 9200A; det ails can be 
found in lYork et alj feOOOl) and IStoughton eFZI j2002l) . Spectro- 
scopic targets chosen by the various SDSS selection algorithms are 
arran ged onto a series of 3" diameter circular fields jSlanton et al.l 



120030. QSO targets are selected on colou r and z— band point sp read 
function (PSF) magnitu de; for details see [Richards et^j2002h and 
ISchneider et alj d2007l) . QSO re dshifts ar e predominantly derived 
from a fit to the QSO template of lvanden B erk et al. ( 200 1 ). 

The third data re lease (DR3) of the sp ectroscopic survey cov- 
ers 37 32 sq. degrees jAbazaiian et al.ll2005) . The DR3 QSO cata- 
logue dSchneider et alj|2005h provides a more complete sample of 
QSOs than the SDSS standard pipeline and corrects the redshifts of 
a small number of objects. 



3 THE ABSORPTION LINE SAMPLES 

The absorption line samples used in this paper are based on 
the absorption line catalogue described by York et al. (2006) 
and York et al. (2008, in preparation). The catalogue contains 
a list of narrow (<700kms~^) absorption features detected 
in QSO spectra draw n from the SDSS DR3 QSO catalogue 
( ISchneideretal]|2005h . Independent searches are carried out for 
both the C IVAA1548.2, 1550.8 and Mg iiAA2796.4, 2803.5 dou- 
blets, and the equivalent width (in the QSO frame) and significance 
of individual line detections are catalogued. We select C IV ab- 
sorbers with equivalent width of the strongest line in the doublet, 
A1548.2, greater than 0.3A. We apply the same equivalent width 
cut to the Mg II doublet. 

In this paper, we retain only those systems with detection sig- 
nificance > 4(7 for the primary line, and > 3(t for the secondary 
line. In each sample, duplicate entries are also removed. These arise 
because the master catalogue contains the result of independent 
searches for Mg II, C IV, and Fe II. We remove duplicates with 



exactly the same Zabs 



Q, and line equivalent widths (EWs). If 



there are two absorbers with redshifts within 0.005 of each other 
(corresponding to AA,- of 7.7A for C IV and 14A for Mg II), the 
strongest system is retained. 

The absorber sample is then cross-matched with the SDSS 



DR3 specobj view in order to obtain a well-defined QSO sam- 
ple with photometric and spectroscopic measurements. The FIRST 
radio survey fluxes are extracted from the DR3 QSO catalogue 
('Bec keretal]|l995l) . This gives a total of 15416 and 24406 C IV 
and Mg II absorption systems. 

Finally, the SNR of the QSO spectra must be sufficient to al- 
low lines above a certain equivalent width to be detected through- 
out their wavelength coverage, otherwise we may introduce arti- 
facts into our correlation functions. Monte Carlo techniques could 
in principle be employed to simulate the full selection function of 
an absorption line search. However, with the large number of ab- 
sorbers and QSOs available to us from the SDSS we can afford to 
make stringent cuts to ensure a (roughly) uniformly selected sam- 
ple. 

We therefore remove absorbers which fall below 4000A or 
above 9000A observed frame wavelength to avoid the noisiest re- 
gions of the SDSS spectra. Figure [T] shows the r-band SNR of the 
QSO spectra vs. EW of the C IV absorption lines (computed in the 
absorber rest-frame). There is a sharp rise in the median and mean 
SNR of the spectrum as the EW of the absorber drops below 1 A: i.e. 
weak absorbers are preferentially found in high SNR spectra. Cut- 
ting at per-pixel-SNR> 8 (where a pixel covers 69 km/s in velocity 
space) and EW> 0.5A gives roughly constant mean and median 
SNR as a function of absorber EW (red and green lines). We there- 
fore adopt these cuts, reducing our sample to 7393 C IV and 16173 
Mg Ilabsorbers. 

We have chosen not to eliminate BALQSOs from our sample, 
as has been done in most previous studies. Our reason for retain- 
ing them is that a BALQSO can also show narrow absorption at a 
different velocity to the broad absorption. As the physical process 
responsible for the presence of a BAL is not well understood, we 
do not wish to bias the NAL sample by removing QSOs with BALs 
from our sample entirely. Unfortunately the absorption line cata- 
logue suffers from a large number of "false positive" line identifica- 
tions in the regions of the BAL troughs. BAL troughs are more of- 
ten observed in C IV than in Mg II so we visually inspected all C IV 
absorber candidates to eliminate false positive detections, leaving 
us with 6456 C IV absorbers. There is a continuum of line widths 
from NALs through to BALs; the criteria imposed during this vi- 
sual inspection were that the lines had velocity widths smaller than 
700 km/s and both members of the doublet were clearly identifi- 
able as an individual system in velocity space, rather than being 
part of broader velocity structures. Except for the most saturated 
cases, the two lines of the doublet are clearly distinguishable using 
this velocity width limit. As we do not determine the complete- 
ness of our absorber sample for every QSO spectrum, the inclu- 
sion of the BALs will cause us to underestimate number densities 
slightly where the flux in the BAL troughs r eaches zero. However , 
ttie fraction o f BALQSOs is around 10-15% faewett & Foltzl2003l ; 
Ijrump et al. 2006) and their typical line widths are one to a few 
thousand km/s. Thus, the fraction of the total pathlength lost due to 
BAL troughs is small. 

For the larger Mg II sample, only absorbers within v/c < 
0.015 of the background QSO were inspected, and only false pos- 
itive identifications caused by BAL troughs and poor continuum 
fits around the QSO emission lines were removed. This procedure 
leaves 16137 Mg II absorbers. In fact, the majority of false positive 
line identifications in the Mg II sample are caused by OH skylines 
(see Wild & Hewett 2005). 

Our samples of C IV and Mg II absorbers have mean redshifts 
of 2.08 and 1.17 respectively. Figure[2]shows the observed redshift 
distribution of the C IV and Mg II samples, before and after the 
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Figure 1. r-band SNR of the background QSO spectrum vs. EW of detected 
C IVA1548 lines. Straight blue lines show the sample cuts used when mea- 
suring the cross-coiTelation of absorbers with QSOs. The dashed red and 
solid green lines show the mean and median SNR as a function of EW for 
QSOs above the horizontal blue line. 



SNR and EW cuts are made. The cuts do not appear to significantly 
bias the redshift distributions. 



4 LINE-OF-SIGHT DISTRIBUTION OF QSO-ABSORBER 
SEPARATIONS 

To get a first impression of the quaUty of our data. Figure [3] shows 
the number density of C IV and Mg II absorption fine systems as 
a function of comoving distance (r) from the background QSO. In 
both cases, the distributions have been corrected for the fact that 
the pathlength available to find an absorber decreases as comoving 
distance increases. The general form of these distributions is ro- 
bust to even significant changes in the sample. Here, we choose to 
present line-of-sight distributions in comoving distance rather than 
the more commonly used velocity parameter /3 = u/c, as we aim to 
quantify the distribution of absorbers caused by galaxy clustering. 
We will return to the {3 distributions in Section[71 

We note that there are clear similarities between the line-of- 
sight distributions of the Mg II and C IV absorber samples. After 
correction for survey completeness, both show a near constant ab- 
sorber density at large distances from the QSO, as we would ex- 
pect for an intervening absorber population. They also both display 
a very pronounced enhancement in the fraction of absorbers within 
25-50 Mpc from the QSO. The C IV absorbers in addition show an 
enhanced number of absorbers out to 150 Mpc from the QSO. A 
similar, although much smaller enhancement may be present out to 
75 Mpc in the Mg II absorbers. 

As will be described more fully in the following sections, the 
observed distributions may be interpreted as the superposition of 
(i) a narrow spike at r < 50Mpc due, in part at least, to galaxy 
clustering; (ii) a tail to high velocities caused by sub-relativistic, 
moderately ionised outflows driven by the QSO (seen primarily in 
the C IV distribution ); (iii) a constant background level of inter- 
vening absorbers. 

Finally, it has been claimed that systematic uncertainties in the 
estimated redshifts of the higher redshift QSOs in our sample could 
potentially lead to false conclusions being drawn about the nature 



of the associated absorbers l lRichards et al.l2002h . The full effect of 
such systematic line shifts on the SDSS redshifts derived by fitting a 
QSO template to the SDSS spectra is complex. A detailed analysis 
of the full redshift-dependent systematic shifts will be presented 
in Hewett & Wild (2008, in preparation). A preliminary analysis 
comparing the line-derived redshifts to the quoted SDSS redshifts 
leads us to shift QSOs with 1.1 < z < 2.1 redwards by 200km/s 
and QSOs with z > 2.1 redwards by 415 km/s. Although shifts of 
this order of magnitude appear warranted by the present dataset, the 
precise values used make only small differences to our final results, 
and are likely to be refined by future work. The dotted histogram in 
Figure|3]shows the line-of-sight NAL distribution before correcting 
the QSO redshifts for this bias, while the solid histogram shows the 
corrected distribution. The difference is very small. 



5 3-D CORRELATION OF ABSORBERS AND QSOS 

In this section we cross-correlate the absorbers with QSOs on 
neighbouring lines-of-sight, in order to estimate the contribution 
to the line-of-sight absorber distribution from galaxies clustered 
around the QSO hosts. 

I t is generally believed that QSO s reside in galaxies 
(e.g. iKristiatJ 1 19731; iDunlop et all Il993h . Strong metal ab- 
sorbers also arise prim a rily i n ordinary galaxy halos and 



disks teahcaU & Spitzeil 1 19691; ISoksenberg & SargentI 1 19781; 



ISteidel et alJll997l ; IZibetti et alj|2007h . Because galaxies are clus- 
tered, we expect to see some enhancement in the number of ab- 
sorbers close to the QSO. With the high density of QSOs available 
in the SDSS survey, one can estimate this enhancement by counting 
the number of absorbers along lines-of-sight that pass close to each 
QSO. 

The standard measures of galaxy clustering are the 2-point 
auto- and cross-correlation functions. These are typically well ap- 
proximated by power laws: 



«(0 = - 
ro 



(1) 



where 7 « 1.8 and rp « 5/i^^Mpc for bright galaxies at low 
redshifts ( |Peebleslll98(lt) . The ratio of the number density of pairs 
separated by distance r to that expected for a random and uniform 
distribution of objects with the same number density is 1 + ^. In 
this section, we estimate ^qa, the QSO-absorber cross-correlation 
function. 0. 



5.1 3-D correlation estimation 

The basic principle of our estimation method is to compute the ex- 
cess number of QSO-absorber pairs relative to an unclustered sam- 
ple, as a function of comoving separation. The number of QSO- 
absorber pairs expected for a uniformly distributed sample is deter- 
mined by computing the number of absorbers expected at random 
along the actual sightlines used to find the observed pairs. In this 

^ We note that in reality we measure the redshift-space cross-correlation 
function. We do not attempt to correct for redshift space distortions: our 
approach is conservative because on small scales, < 10 Mpc, non-linear z- 
space distortions will act to reduce the observed clustering amplitude in a 
similar way to the effect of eiTors on the QSO redshifts. On larger scales, 
the effect of linear z-space distortions o n the correlation fu nction is small 
at the redshifts of interest here (see e.g. ICroom et al.ll200 for a detailed 
discussion). 
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Figure 2. Left: The redshift distribution of the C IV absorbers before (full) and after (dotted) the cuts on QSO spectrum SNR and absorber EW. Right: The 
same for the Mg II sample. 
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Figure 3. Left: The number density of line-of-sight C IV absorbers found in QSO spectra as a function of comoving distance from the QSO, after (full) and 
before (dotted) correction of the QSO redshifts for the bias caused by blue-shifted UV emission lines. The distribution has been corrected for completeness 
due to the decreasing pathlengths available to find an absorber as distance increases. Note the peaked excess of absorbers lying very close to the QSO, the 
significant number of negative velocity absorbers, and tail of excess absorbers extending to inferred distances of around 150 Mpc. Right: The same for an 
independent sample of Mg II absorbers. An excess is again observed close to the QSO; a small tail to larger distances may also be present. 



way, we automatically account for selection biases such as those 
caused by fibre collisions in the SDSS survey. Our final measure- 
ment yields the excess number of absorbers as a function of co- 
moving distance from a QSO. We can then use this measurement 
to ascertain whether the excess seen along the line-of-sight is con- 
sistent with galaxy clustering alone, or whether some component 
of the excess is intrinsic to the QSO itself. The scheme is presented 
pictorially in Figure|4] to which we will refer throughout this Sec- 
tion. 

First we build our QSO-QSO pair catalogue. For each QSO 
in the DR3 QSO catalogue, the SDSS Catalogue Archive Searcher 
(CAS) is used to find all QSOs at higher redshift and within 3 de- 
grees of it on the sky. We cross-match the DR3 QSO catalogue with 
the specobj view of the CAS in order to ensure a unique sample of 
QSOs from main SDSS survey plates. The DR3 QSO catalogue 
contains 45,988 QSOs; this results in ^^6.5 million pairs above a 
redshift of 0.3. Only 10823 DR3 QSOs have SNR above the thresh- 



old (8 per pixel), and Zqso > 1.58 so that C IV can be seen above a 
restframe wavelength of 4000A. This reduces the number of QSO- 
QSO pairs drastically for the C IV analysis compared to the Mg II 
analysis. QSOs below the SNR threshold are retained only as cen- 
tral QSOs (see Figure|4ll in the correlation analysis, i.e., the sight- 
lines provided by their spectra are not included in the absorption 
system search or in the control absorber counts. 

Next, we select the sightlines which we will use in the analy- 
sis. We only consider sightlines from QSOs above the SNR thresh- 
old. As indicated in Figure |4] we restrict the sightlines to start at 

> 0.1 from their QSO host, and to lie redwards of 1250A in 
the QSO rest frame (to avoid the strong Lyman-a line at 1215A 
and the Lyman-a forest bluewards of Lyman-a). The velocity dif- 
ference criterion Az > 0.1 ensures that the final signal is not af- 
fected by the excess of associated absorbers with sub-relativistic 
velocities described in Section |4] Because of the restriction that 
sightlines must end at 1250A we cannot increase this limit much 
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Figure 4. A pictorial representation of the method used to calculate the 3-D cross correlation between QSOs and absorbers. From a central QSO we count, 
using small bins of A^, the observed and expected number of absorbers along neighbouring sightlines to background QSOs. These numbers are accumulated 
as a function of the comoving separation, r, between the center of the Az bin and the central QSO. The process is repeated for all central QSOs, thereby 
building a 3-D cross-correlation estimate. 



further without drastically reducing the number of QSO- C IV ab- 
sorber pairs. Contained within these sightlines are 3374 C IV and 
13504 Mg II absorbers. 

Using these sightlines and absorbers, we calculate the number 
density of absorbers as a function of redshift n(z) shown in Fig- 
ure [5] and use this to estimate the expected number of absorbers in 
each bin along each sightline. Then, for each central QSO we 
use the sightlines of all its paired QSOs to measure the expected 
and observed number of absorbers as a function of comoving sep- 
aration between the central QSO and the position along the sight- 
linfl In Figure|4]the arrow, labelled r, indicates the comoving dis- 
tance between the central QSO and one of the bins along one of the 
sightlines. In this bin, no absorber is found so no addition is made 
to the observed number of absorbers at this comoving separation. 
However, a small addition will be made to the expected number of 
absorbers at this separation, based on the pre-calculated n{z) at the 
redshift of this Az bin. As well as accruing the observed and ex- 
pected numbers of absorbers as a function of comoving separation, 
we sum the number of contributing sight lines, which is necessary 
in the estimation of the errors. Finally, because we must bin our 
results in Ar, and we do not have enough sightlines to create in- 
finitesimal bins, we must account for the fact that more sightlines 
within a Ar range will lie towards the outer edge of the range than 
the inner edge i.e. the effective radius of our Ar bin is not centered 
on to the bin, but biased towards the outer edge. We do this by ac- 
cumulating the expected number of absorbers, weighted by r~^'^ 
(see Equation|3]l. 

The final 3-D correlation estimate is given by: 



1 + e(rcff ) = 



No(r,ff) 
No(roff) 



(2) 



^ In a flat universe, the comoving distance between two objects is given by 
the square root of the sum of the squares of the angular and line-of-sight 
comoving distances. 



where No(rcff) is the observed number of absorbers in a Ar bin 
with effective radius rcjr, and Nc(roff) is the expected number as 
described above. The effective radius of each of the Ar bins is 
estimated to be: 

'j^Ner-^x{-V7) 



Toff 



(3) 



where the sum is over all Az bins (segments of the sightlines in 
Figure|4]l that contribute to the Ar bin. 7 is chosen to be 1.7, which 
is close to the measured exponent in the Mg II-QSO correlation 
(see below), although the precise value used does not change the 
results significantly. 

On small scales the QSO-absorber pairs in each Ar bin are 
independent (in general, absorbers are only paired to one QSO on 
small scales), thus the error is estimated assuming Poisson statis- 
tics: 



Ae(r) 



1 + CW _ x/TO 



Nc(r) 



(4) 



On larger scales the pairs in each Ar bin are no longer indepen- 
dent (one absorber can be paired with many QSOs in the same 
bin ) and the Poisson errors underestimate the true errors. Follow- 
ing [shanks|&|Boy3 il994h . as the number of pairs approaches the 
number of absorbers in the analysis (Nabs), the errors are approxi- 
mated by: 



Ae(r) = 



i + e(r) 

N/(Nabs) 



(5) 



Finally, we calculate the significance of a positive clustering 
signal (excess number of pairs) assuming a Poisson distribution. 
Assuming no clustering, the probability of observing No or more 
pairs at any given r is given by: 

P(fc>N,)= Y. ■ ^6) 



fc! 
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In the closest bins with fewest observed pairs, we checked that the 
probabilties derived here agree with those derived from a Binomial 
distribution. 



5.2 Clustering detection 

The correlation functions for C IV and Mg II absorber-QSO pairs 
are shown in Figure [6] where we plot (1 + in bins of ~ 6 Mpc. 
The results are plotted logarithmically in the right hand panels. As 
expected, there is an enhancement in the counts of both QSO-C IV 
pairs and QSO-Mg II pairs at small comoving distances from the 
QSO. The significance of the detection is considerably higher for 
the Mg II absorbers, because the sample is much larger. Overplotted 
as a dashed line is the 68% confidence limit for detection of a clus- 
tering signal, given the observed and expected number of absorbers 
in each bin (Eq.[6](. In other words, if there were no clustering, we 
would expect an independent bin to lie above this line 32% of the 
time. 

We can evaluate the amplitude and significance of the clus- 
tering signal for the C IV systems as follows. The total fractional 
excesflfl of C IV absorber-QSO pairs between ~6 and ~43 Mpc 
(bins 2-6) is 1.49. From Eq.[6l the probability of obtaining a value 
greater than or equal to this value is 0.006 i.e. we detect excess 
clustering at a 99.4% confidence level. If we instead use only the 
second to fifth bins, we obtain a similar confidence level of 99.7%. 
Including the first bin results in a similar significance detection, but 
a lower clustering amplitude. We have chosen to discard this first 
bin, which only contains one pair, because of the uncertainty in 
the QSO redshifts for this sample. Redshift errors will act to move 
clustered QSO-absorber pairs to outer bins, and cause a flattening 
in the correlation function at small distances. Our choice to discard 
the first bin results in a small overestimation of the clustering am- 
plitude, as some QSO-absorber pairs from the central bin will be 
counted in outer bins. The number of QSO-absorber pairs expected 
at these small radii is very low, however, and their contribution to 
the derived clustering amplitude is not significant within the final 
errors. 

For the Mg II absorber sample, the fractional excess within 
43 Mpc is 4.0 and the detection of clustering is significant at the 
>99.9999% level. 

These fractional excess numbers of QSO-absorber pairs (1.49 
for C IV and 4.0 for Mg II) can be fit to the standard power law form 
used to describe the correlation function, ^ (Equation^. To be con- 
sistent with typical values quoted in the literature, we calculate ro 
forafixed7 = 1.8. We find ro = 5.8±1.1 and 4.84±0.4/!,"^ Mpc 
for the C IV and Mg II absorbers respectively. These fitted power 
laws are overplotted in Fig.|6] For the Mg II absorbers we can also 
fit the correlation function to all bins directly, leaving 7 as a free 
parameter. We use a non-linear least squares technique to obtain 
7 = 1.67±0.09 and ro = 5.0±0.4/i"^ Mpc. The resulting model 
has a reduced of 1 .007. The curve is overplotted as a blue dotted 
line in Fig.|6] In conclusion, in this section we have unambiguously 
detected and measured the excess of absorbers in the vicinity of 
QSOs due to the clustering of galaxies. 



6 A MODEL FOR LINE-OF-SIGHT NALS 

In this section, we model the line-of-sight correlation of QSO- 
absorber pairs and we address two questions: (i) is the central spike 
at Zabs ^ ZtiBo consistent with being caused by galaxy cluster- 
ing? (ii) can the sub-relativistic tail of C IV absorbers be caused 
by galaxy clustering? 

We proceed by constructing a model with three components: 
(i) the amplitude of the clustering of absorption line systems arising 
from galaxies in the vicinity of the QSOs; (ii) the velocity offset and 
dispersion induced by peculiar motions and redshift errors; (iii) the 
distance to which the QSO ionises surrounding gas. 

Motivated by the results presented in the previous section, we 
construct a simple toy model of the form: 



nabs{r) = 



1 + 







r ^ -Rcut 



■■g{r,(j) (7) 



^ Fractional excess is defined as (Nci-Nc)/Ne = 
© 0000 RAS, MNRAS 000, 000-000 



where Uabs is the mean number of absorbers per comoving Mpc at 
comoving distance r from the QSO, rihg is the comoving number 
density of intervening absorbers, -Rcut is the comoving radius out 
to which the QSO photoevaporates the absorbers (the "proximity" 
zone), and ro is the correlation length of QSO-absorber clustering. 
The model is convolved with a Gaussian of width a to reproduce 
the effects of the departure of the QSOs and absorbers from the 
Hubble flow, as well as errors on the redshift measurements of the 
QSOs. 

Of the four parameters, ro is measured from the 3-D correla- 
tion estimate (Section [5]l and n^g is measured at distances greater 
than 200 Mpc from the QSOs. There are two remaining parameters: 
a and -Rcut . We will discuss plausible values for these parameters 
in the following subsections. Our model is no doubt a simplifica- 
tion of the true situation. For example, we do not explictly account 
for a contribution from the halo clouds around the QSO, assuming 
only a continuation of the powerlaw component towards the QSO. 
In SectionQwe will discuss the limitations of the model. 



6.1 Velocity dispersion 

There are two contributions to the dispersion o: redshift errors on 
the QSO and, assuming the QSOs and absorbers originate in dif- 
ferent galaxies, peculiar motions of the absorbers and QSO host 
galaxies with respect to the Hubble flow. 

The SDSS catalogue provides errors on the redshifts of each 
QSO. The median quoted error of all the QSOs with absorption line 
systems are 0.0017 for the Mg II absorber sample and 0.0019 for 
the C IV absorber sample. This corresponds to velocity errors of 
180km/s for both samples at the median QSO redshifts (1.89 and 
2.19 for the Mg II and C IV samples). However, given the observed 
syst ematic shifting of broad QSO lines from the systemic redshift 
(e.g. lGaskell|[T98a ; ISulentic et ai]|200a) . the true errors are likely 
to be considerably lar ger, particularly for th e higher redshift C IV 
sample. As shown bv iRichards et alj ( l2002 l) C IV is shifted by a 
median ~800km/s from the lower ionisation Mg II line. The shift 
can be as much as 3000 km/s in some QSOs. The position of the 
C IV line is used in the template to derive the systemic redshift of 
the SDSS QSOs and we may therefore expect this shifting of the 
lines to cause errors of up to ~ 1000 km/s. 

The second contribution to the dispersion o comes from the 
peculiar motions of galaxies with respect to the Hubble flow, if the 
absorber and the QSO reside in different hosts. The pairwise veloc- 
ity dispersion of galaxies ((T12) has been measured at low redshift 
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Figure 5. Left: Number density of detected C IV absorbers as a function of Zabs> niz). Overplotted in red is the smoothed distribution used in the analysis. 
Right: Same for Mg II absorbers. 
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Figure 6. Top: The 3-D correlation between C IV absorbers and QSOs (i.e. observed over expected number of pairs), as a function of comoving distance 
from the QSO. The right-hand panel is a logarithmic rendering of the left-hand panel. The overplotted plain line is the powerlaw described in the text, with 
fitted ro = 5.8 it l.lh~^ Mpc for fixed 7 = 1.8. The dashed line is the 68% detection threshold given the number of sightlines and number density of 
absorbers {Eq.[6). Bottom: Same for Mg II absorber-QSO pairs. The additional dotted line is the fully fitted powerlaw, with ro = 5.0 ± 0.4/i~^ Mpc and 
7 = 1.69 ± 0.09. 
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in th e 2dFGRS and SDSS surveys l ljing & B6merll2004l ; llTetai] 
120061) . It can range from a few h undred to near ly a thousand km/s, 
depending on the type of galaxy jLi et al .120061) . However, the red- 
shift evolution of ai2 has not been measured, making it difficult to 
predict the velocity dispersion experienced by galaxies in the rela- 
tively massive haloes of the QSOs. 

As neither the redshift errors nor the expected peculiar veloci- 
ties of the galaxies are well known, we choose to estimate a directly 
from the data, by ensuring that our models reproduce the observed 
width of the central peak in the comoving separation distribution 
(Figure [3]l. The derived scatter in comoving distance is 9.0 Mpc 
for the Mg II absorber sample and 12.2 Mpc for the C IV absorber 
sample. This corresponds to velocity dispersions of 580 km/s and 
1030 km/s at the median redshifts of the Mg II and C IV absorbers 
in the central peak of the comoving separation distribution (1.35 
and 1.97, see Section [TJ. 

It is important to note that the precise value chosen for a does 
not greatly affect our results; it only slightly alters the shape of the 
central peak in the distribution of absorbers in velocity space. 

6.2 The line-of -sight proximity effect 

Two further physical scales are important for the geometry of line- 
of-sight QSO-absorber clustering. The first is the distance to which 
the gaseous halo of the QSO itself can contribute to the number 
of pairs. The second is the distance out to which the intense radia- 
tion field of the QSO eliminates the particular ion in question. Here 
we simply present some typical scales culled from the literature 
to guide in the interpretation of our results, J?cut is left as a free 
parameter in our model. 

As discussed in the introduction, halos of Mg II around galax- 
ies have been studied for many years. Their proper e xtent is 
believ ed to be ~40/i~^kpc for unity covering fraction jSteidell 
Il993h . Beyond this distance absorbing clouds have been found 
out to proper radii of ~ 100 kpc, but wit h covering factors 
that decrease with dis tance CBowen et al. 200d: Izibetti et alJl 20071 : 
iKacprzak et alj|2007h. Fewer q uantitative studies have been made 
for C IV. but lChen et al.l ^200 11) find C IV halos of ~ lOO/i"^ kpc 
in proper radius, with unity covering factor. C IV and Mg II are 
believed to originate in the same halo clouds, with the higher ion- 
isation ion C IV alm ost always found in conjunction with Mg II, 
but not the converse ^Bergeron & Boisse|[l984l : ISteidel & Sargent! 
ll992l : IChurchill et aljl999r The true extent of absorber clouds de- 
pends sensitively on the equival ent width limit imposed, also possi- 
bly on the mass of the host halo jSteidel et aljl994l : lchen & TinkeJ 
l2008i) and it is likely that it changes with redshift. We therefore take 
these nominal values for convenience, but the reader should keep 
in mind that they are only indicative. 

The effect of the ionising field of the QSO on the surround- 
ing inter -galactic medium is less well constrained ob servationally 
(see e.g. lEUingsonet afll 1994b . IChelouche et"^ J20"07i) have used 
observational constraints from Mg II absorbers to model the effect 
using the ionisation code CLOUDY jFerland et al.lll99^ . Assum- 
ing a uniform radiation field and a cloud density n = 0.02cm~"^, 
they conclude that Mg II will be completely destroyed in the in- 
ner 100 kpc of the halo and that the density of Mg II absorbers 
will only recover to the background level on Mpc scales. They also 
predict that the fraction of C IV systems will fall off rapidly be- 
low ~200kpc. If the QSO radiation field is anisotropic, the line-of- 
sight ionisation boun dary must move further out to remain consis- 
tent with observation. iHennawi & Prochaskal j2007h also present a 
simple model of the proximity effect for Lyman-limit systems. As 



Mg II systems are found to correspond closely to Lyrnan-lim it sys- 
tems teergeron & StasinskalfTgsa : ISteidel & Sargenllll992h . these 
results are relevant to our Mg II sample. They find absorbers with 
neutral hydrogen column densities below lO^^atoms/cm^ and hy- 



drogen volume densities below 0.1cm 
within 1 Mpc of the QSO. 



to be photo-evaporated 



7 RESULTS 

In Table [T] we present the fraction of the QSOs in our sample that 
have absorbers with —0.01 < /3 < 0.01 (associated) and with 
0.01 < 13 < 0.04 (distant). 3.4% (15%) of QSOs in our sam- 
ples have one or more associated Mg II (C IV) absorbers with rest 
equivalent width > 0.5A. We note that these samples are at dif- 
ferent mean redshifts, and th erefore the compari son of the fractions 
must be treated with caution. IVestergaarj ( |2003|) find 25% of QSOs 
to have C IV with the same EW limit within ±5000 km/s. We find 
18% within this higher velocity limit, which is con sistent given the 
small number of objects in the I Vestergaarc) ( l2003h sample. Given 
the fact that Mg II and C IV are almost always seen if a QSO sight- 
line passes close to a galaxy, these results indicate that the QSO 
must ionise most of the Mg II and C IV clouds in its own halo. 

In this section it is often necessary to convert between physical 
and comoving scales. As we are primarily interested in describing 
the central peak of the comoving separation distribution, we use the 
median redshift of those QSOs contributing to the peak. These are 
1.35 and 1.97 for the Mg II and C IV samples respectively. 



7.1 The Az ~ excess 

We now compare the line-of-sight distributions of the NALs with 
the prediction of our model for the component due to galaxy clus- 
tering. We first assume that the QSO ionises clouds only in its own 
halo, i.e. we set iicut ~ Rhaio in Equation [7] We take i?haio to 
be the comoving size of absorber halos at the median redshifts 
given above assuming typical sizes from the literature, i.e. 130 and 
420 kpc for the Mg II and C IV samples respectively. The predicted 
line-of-sight absorber distributions are shown in Figure|7] There are 
two clear results, each with an obvious implication. Firstly, galaxy 
clustering only contributes to the central ~ 50 Mpc of the distribu- 
tion. We therefore conclude that the tail of high- velocity absorbers 
seen in the C IV distribution must be intrinsic to the QSO and/or 
its host galaxy. Secondly, the model greatly overpredicts the num- 
ber of Mg II absorbers, leading us to conclude that the QSOs ionise 
Mg II clouds well beyond the halo of the host galaxy. The model 
overpredicts the Mg II counts within ±40 Mpc of the QSO (corre- 
sponding to ±2500 km/s) by a factor of 2.05. 

By increasing the size of the region over which the QSO 
ionises Mg II (Rcut), our simple model can reproduce the ob- 
served distribution of Mg II absorbers reasonably well (albeit with 
a small offset that we will address below). Figure[8]shows the line- 
of-sight distribution of Mg II and C IV absorbers within 200 Mpc 
from the QSO. Overplotted in the right hand panel is a model with 
Rcut = 0.8 Mpc (close to the size of the proximity zone for Lyman- 
limit systems derived by Hennawi & Prochaska [2007]). 

Although our model with a large proximity zone appears to fit 
the Mg II data rather well, we note that our simple "top-hat" exclu- 
sion zone is no doubt an over-simplification. Some Mg II absorbing 
clouds probably do exist within 0.8 Mpc of the QSO. Such systems 
may be associated with material within QSO host galaxies. Vanden 
Berk et al (2008) have shown that such associated Mg II systems 
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Table 1. The fraction of QSOs with one or more strong Mg II or C IV absorption lines within fixed velocity ranges. The absorber rest frame EQW limit is 0.5A 
and the radio-quiet (RQ) and radio-loud (RL) samples have been matched in K-corrected i-band absolute magnitude. 



velocity FMgll(all) FMgll(RQ) FMgll(RL) FciV{all) FciV{RQ) Fciv{RL) 



-0.01 < /3 < 0.01 


0.034 


0.036 


0.054 


0.15 


0.14 


0.16 


0.01 < 13 < 0.04 


0.044 


0.051 


0.055 


0.11 


0.11 


0.12 
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Figure 7. Left: The number of C IV absorbers per comoving Mpc as a function of line-of-sight comoving distance from the QSO. Overplotted as a dashed line 
is the modelled galaxy clustering component, assuming the QSO photoevaporates a sphere corresponding to previous estimates of the typical size of absorber 
halos. Right: same for Mg II absorbers. 



tend to be on average dustier (implying denser gas) and be more 
highly ionised (as would be expected for clouds within Rcut)- By 
increasing the radius out to which Mg II halo clouds of average 
volume density are evaporated, our model would allow for an in- 
creasing fraction of the associated absorbers to be directly linked to 
the presence of the QSO and its host galaxy. We will return to this 
point in Section[T2] 

Unlike the Mg II systems, the C IV absorber distribution does 
not require a large proximity zone. Figure[7]shows that if the QSO 
photoevaporates all C IV clouds in its own halo, the number of ab- 
sorbers within ±40 Mpc is underpredicted by more than a factor of 
two. We can try to fit the distribution by decreasing Rcut ; the model 
over-plotted as a dotted line in the left-hand panel of Figure [8] has 
Rcut =0.18 Mpc, corresponding to a proper radius of 60 kpc at the 
median redshift of the QSOs. This radius is smaller than the esti- 
mated, but uncertain, C IV halo size discussed in Section|4] We note 
that such a small proximity radius may conflict with the small frac- 
tion of QSOs observed to have associated C IV absorption (16%). 
We will return to this point in the following sections. 

The main conclusions of this section as follows: 



(i) the central peak of the Mg II absorber distribution could be 
explained by galaxy clustering if the QSO ionises absorbers to ~ 
0.8 Mpc. Assuming the physical size of Mg II halos remains the 
same as observed at low redshift, this corresponds to more than 6 
times the size of the QSOs own Mg II halo, implying that absorbers 
in the halos of neighbouring galaxies along the line-of-sight are 
also affected. 

(ii) The ionising sphere must be substantially larger for the 
Mg II absorption line systems than for the C IV absorption line 



systems, consistent with the higher ionisation energy of the latter 
ion. 

(iii) The C IV systems exhibit a high velocity tail of absorbers, 
which cannot be explained by galaxy clustering and hence must be 
intrinsic. It is likely that these intrinsic absorbers also contribute to 
the distribution at lower velocities. 

In the following subsection we will further investigate the dis- 
tribution of these intrinsic absorbers. 

7.2 The high velocity C iv absorbers 

From Figure [8] it is clear that many narrow C IV absorbers are in 
a high velocity tail that extends out to 150 Mpc before becoming 
indistinguishable from the background absorber population. This 
tail cannot be due to galaxy clustering, and we conclude that these 
absorbers are truly intrinsic to the QSO. 

In Figure [9] we present a simple model for the velocity dis- 
tribution of the absorbers. We have fit an exponential to the high- 
velocity component: 

Wabs - ntg cx e""'^ (8) 

where 

c R'' + 1 
R^\±^, (10) 

1 + Zabs 

and a = 70 ± 7.7. The model is fit in the range 0.01 < f3 < 0.05 
to avoid the central clustering peak. 

Our model predicts that between 3000 and 12000 km/s 
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Figure 8. Left: The number of C IV absorbers per comoving Mpc as a function of line-of-sight comoving distance from the QSO, focussing on the region 
close to the QSO. Overplotted as a dashed (blue) line is the modelled galaxy clustering component, assuming the QSO photoevaporates a sphere of radius 
O.lSMpc, smaller than the observed size of C IV halos at low redshift. The full line is for the full model with Rent = 0.3Mpc and including the intrinsic, 
outflowing absorber distribution. The dotted (red) line is the clustering component of this final model. Right: same for Mg II absorbers, with the clustering 
model assuming a proximity sphere of radius 0.8 Mpc, more than six times larger than the Mg II haloes observed at low redshift, and assuming no intrinsic 
component. 



(0.01 < /3 < 0.04), 45% of C IV absorbers are caused by sub- 
relativistic outflowing material from the QSO. Due to the presence 
of the clustering component, we cannot constrain the shape of this 
intrinsic component at low- velocities. While it seems unlikely that 
the number of intrinsic absorbers drops sharply towards low ve- 
locities, various processes may act to enhance their number (see 
Section rot . Our derived fraction of 45% intrinsic absorbers may 
thus be an underestimate. 

Along similar lines, we note that a small excess of Mg II ab- 
sorbers with 0.005 < (3 < 0.015 may be present in Figure[8] We 
have attempted to constrain the strength of this excess by fitting an 
exponential distribution with a width fixed to be that observed in 
the C IV absorbers. We obtain a height which is 20±5% that found 
for the C IV absorbers. 

We now attempt to build a model to describe both the extrinsic 
and intrinsic C IV systems. We convert our modelled velocity dis- 
tribution of outflowing material into a distribution in comoving dis- 
tance by assuming that each QSO contributes an equal proportion 
of the high-velocity tail. For each absorber-QSO pair we transform 
the model velocity distribution into the distribution in comoving 
distance at the redshift of the QSO. We convolve the final distri- 
bution with a Gaussian to account for the QSO redshift errors and 
add this to the clustering model with different values of the param- 
eter Rcut- The result with iicut = 0.3 is shown in Figure[8]as the 
solid black line. The contribution from clustering is shown by the 
dashed red line. The model fits well and it implies that a substan- 
tial fraction of low-velocity C IV absorbers, with v < 3000 km/s, 
may also be part of a distribution of outflowing material that is in- 
trinsic to the QSO and/or its host. This increased value for iicut is 
now about 70% that of low redshift C IV halos. Together with the 
likelihood that Ticut varies with redshift, QSO luminosity and halo 
mass, the fact that only 16% of QSOs show associated absorption 
does not appear incompatible with our model. 

In Table [2] we present the fraction of each component (back- 
ground, galaxy clustering, intrinsic) present in our model in differ- 
ent velocity ranges for both Mg II and C IV absorbers. The number 
of absorbers associated with intervening galaxies in our model is 



the sum of the background and galaxy clustering numbers. We note 
that the fractions of intrinsic absorbers presented here may be lower 
limits; we will return to discuss this in Section[8] 

7.3 Radio Loud and Radio Quiet QSOs 

A key point of interest is whether the line-of-sight distribution of 
absorbers differs between radio-loud (RL) and radio-quiet (RQ) 
QSOs. The S PSS QSO catalogue has been matched with the FIRST 
radio survey teeckeret all 19951). In this section we have included 
only those QSOs in the FIRST survey footprint. RL QSOs have 
been defin ed to be those w ith 1.4GHz luminosities greater than 
lO^'^W/Hz jMiller et aljT99 0). The FIRST radio survey is flux lim- 
ited to 1 mjy, which corresponds to ^ lO^^W/Hz at a redshift of 2, 
close to the mean redshift of the C IV absorbers. We therefore ex- 
pect our RL sample to be complete in the Mg II sample, but some 
RL QSOs in the C IV absorber sample will contaminate the RQ 
sample. As the fraction of RL QSOs is low, around 10% of all 
QSOs, this will not significantly bias the RQ distribution. In the 
Mg II sample we have 1298 RL and 13179 RQ QSOs. In the C IV 
sample we find 543 RL and 5355 RQ QSOs. 

It is know n that a correlation exists between optical and radio 
luminosity (e.g. ICirasuolo et alj2003l : ISikora et al.ll200 7'). We have 
therefore matched the optical luminosities of our RQ QSO sample 
to those of our RL QSO sample to avoid any biases. Optical lumi- 
nosity is defined from the i-band PSF magnitude, K-corrected to 
2 = 2 following [Richards et alj ilOOdt) . This matching procedure 
reduces the RQ samples to 8343 and 3528 for the C IV and Mg II 
samples respectively. The velocity distributions are shown in Fig- 
urefTOl 

Comparing the RL and RQ samples, we see clear differences. 
A narrower and more pronounced spike at /3 ~ is seen for the 
RL QSOs for both C IV and Mg II absorbers. This result for C IV 
has been presented for many fewer QSOs by , Ves tergaard 
We also note that the subrelativistic tail of C IV absorbers is visible 
in both the RL and RQ samples, but is clearly more pronounced in 
the RQ QSOs. This high-velocity tail was observed previously in 
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Figure 9. The /3 = d/c distribution for C IV and Mg II. Overplotted as a daslied line on the C IV high- velocity tail is an exponential curve e with a = 70. 
A curve of the same width is fitted to the Mg II high- velocity tail, with a height that is 20±5% that observed in the C IV absorbers. 



Table 2. The fraction of Mg II and C IV absorbers above the absorber rest frame EQW limit of 0.5 A attributed to background (bg), galaxy clustering (gc) and 
intrinsic (int) components in our final models. A comoving distance of 40Mpc corresponds approximately to a /3 of 0.008 (0.014) at the median redshift of the 
QSOs in the Mg II (C IV) central peak. A comoving distance of 170 Mpc corresponds approximately to a /3 of 0.037 (0.07). 



velocity 


FMgll(bg) 


FMgll(gc) 


FMgll(int) 


FciV(bg) 


FciV(gc) FciV(int) 


-40 < r/Mpc < 40 


0.46 


0.54 


0.0 


0.13 


0.48 0.39 


40 < r/Mpc < 170 


0.99 


0.01 


0.0 


0.53 


0.01 0.45 



RQ QSOs by lRichards et alj ( 1 19991) and lRichardj ( 1200 ih . We will 
discuss these results in more detail in Section [831 



8 DISCUSSION 

Associated absorption line systems account for a small fraction of 
all metal absorbers detected in spectroscopic QSO surveys; how- 
ever, their study is important for understanding the effect of QSOs 
on their host galaxies and their local environment. In the follow- 
ing subsections we will address the implications of our results for 
understanding the nature of associated narrow absorbers and for 
absorption line systems in general. 

8.1 Intrinsic outflowing gas 

The primary aim of our paper is to quantify whether some narrow 
absorption line systems are intrinsic to the QSO and its host galaxy. 
From our analysis alone, there is clear evidence that a high fraction 
of C IV absorbers are intrinsic to the QSO. 45% of C IV systems 
with velocities in the range 0.01 < [3 < 0.04 are intrinsic and are 
well described by an exponential velocity distribution. A similar, 
but considerably smaller fraction of high-velocity Mg II absorbers 
are attributable to the QSO. The high-velocity Mg II tail is 20±5% 
the height of that observed in C IV. Starburst-driven winds do not 
reach such high velocities; it is therefore clear that these absorption 
line systems are a direct cons equence of a QSO in duced outflow 
(see compilation in Figure 2 of lTremonti et al.l2007l) . 

The physical radius at which the intrinsic C IV NALs originate 
cannot be constrained from these results alone. It is possible that 
they originate from the wind close to the accretion disk of the cen- 
tral AGN itself, although it may be difficult to understand the nar- 



rowness of the lines under these circumstances. A key constraint is 
that narrow intrinsic absorbers are observed in both the RL and RQ 
QSOs, contrary to BALs, t he incidenc of whi c h decreases sharply 
as radio power incr eases ( IStockeetal.lll992l ; iBecker et alj|l997l; 
IShankaretal]|2008l) . One possible explanation is that the strong 
X-ray flux observed in RL QSOs limits the self-sh ielding of the 
wind s and thus the widths of the absorption features ( IMurrav et al.l 
In this scenario NALs are simply failed BALs. Alternatively, 
NALs may be due to the viewing angle of the accretion driven wind 
l lElvisl ^O); the high velocities observed in this work will allow 
tight constraints to be placed on the wind geometry. 

It is also possible that some of these absorbers, particularly 
those with lower velocities, may not be directly associated with 
the accretion disk winds, but originate further out in the ISM 
of the galaxy. The QSO spectra contain significantly more infor- 
mation than we have used in the current paper: for example, it 
is possible to derive the ir dust content and ionisation states (see 
IVanden Berk et al J 12008*) and also to look for fine structure lines 
and line variability in multi-epoch observations. Stacking the spec- 
tra as a function of velocity distance from the QSO will increase 
the signal-to-nois e of the weaker lines needed to o btain accurate 
column densities jWild et alj200^ ; lYork et al.l2006h . All of this in- 
formation should enable us to constrain more precisely the physical 
origin of the C IV outflows, just as it has allowed Vanden B erk et al.l 
(2008) to suggest an origin for some associated Mg II systems in 
the gas of the QSO host galaxy. 



8.2 The proximity effect of the QSO 

Our results show that the QSO has a considerable impact on NALs 
along the line-of-sight, ionising ordinary Mg II halo clouds well be- 
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Figure 10. Left: the fraction of C IV absorbers relative to the background level found in RQ (full line) and RL (dashed line) QSOs as a function of velocity 
shift from the QSO. The peaked excess close to the QSO is more pronounced in the RL QSOs but the relativistic tail of absorbers exists in both RQ and RL 
QSOs. The long-dashed line is positioned at unity for reference. The total number of QSOs in each sample is given in the top right. Right: Same for Mg II 
absorbers. 



yond the expected Mg II halo size of the QSO host galaxy. The QSO 
also destroys a considerable fraction of the C IV clouds within the 
host halo. Our model based on galaxy clustering overpredicts the 
number of Mg II absorbers within 40 Mpc of the QSO by a factor 
of two. This is somewhat less than the factor of 4-20 estimated by 
iHermawi & Prochaskal ( l2007h . who studied the distribution of 17 
Lyman-limit systems around QSOs. 

The effective radius of the first bin in our 3 dimensional Mg II 
absorber-QSO cross- correlation analysis is 4.2 Mpc. Our results 
can thus t est whether there is a "t ransverse proximity effect" on this 
scale (e.g. lOoncalves et alj2007l) . We fit the Mg II correlation func- 
tion (Equation [l), this time excluding the central bin. The result- 
ing power-law predicts 38 pairs in the central bin, compared to the 
34 we actually observe. Thus, we find no evidence for a deficit of 
Mg II absorbe r-QSO pairs at this r a dius. This is in agreement with 
the results of Bowen et alj ( |2006|) , iHennawi" & Prochaskj feOOTh 
and jTvtler etal .l2007h . who find no transverse proximity effect for 
Mg II, Lyman Limit systems, or for metal-line systems in general. 



8.3 Radio loud vs. Radio quiet 

Unfortunately our sample of RL QSOs is too small to allow us to 
measure the 3D correlation function of RLQSO-absorber pairs, so 
we can not assess the origin of the larger excess of low-velocity 
absorbers in RL QSOs, visible in both the Mg II and C IV distribu- 
tions. 

There are several possible explanations for the excess number 
of low-velocity absorbers seen in the RL QSOs: 

• the radio jets themselves drive low-velocity, lower ionisation 
or clumpier outflows, 

• RLQSOs ionise their surrounding medium to smaller dis- 
tances than RQQSOs, 

• the RLQSOs are in denser environments than RQQSOs. 

At low redshift, evidence points to wards R LQSOs existingin 
denser environments than RQQSOs (e.g. lSmith & Heckman 1 990h . 
although these low-z RLQSOs have lower luminosities than those 
in our sample. When compared to galaxies of the same mass, RLQ- 
SOs are preferentially found in galaxies at the center of groups and 



halos teest et alj|2007h . iKauffmann et alj ( |2007|) find that the lo- 
cal density of galaxies around RL-AGN is a factor of two higher 
than around RQ-AGN matched in velocity dispersion, redshift and 
stellar mass. All these results, combined with our excess of low- 
velocity absorbers in the RL QSOs, might suggest that RL QSOs 
exist in denser environments. 

On the other hand, radio jets have be en found to drive low - 
velocity (~100km/s) outflows of hydrogen jMorganti et alj|2007h . 
A similar conclusion has been reached for high-redshift radio 
galaxies by Nesvadba et al. (2007). Hydrodynamical simulations 
show how cool gas which may lead to absorption li ne systems ca n 
build up behind the bow shock of radio jets (e.g. lKrausell2002h . 
Such a scenario should be testable, as it would imply a clear rela- 
tion between the presence of absorbers and the orientation of the 
jet. 



8.4 Absorber and QSO halo masses 

Our derived values of the cross-correlation lengths of QSOs with 
C IV and Mg II systems (ro = 5.8 ± 1.1 and ro = 5.0 ± 
0.4/t^^ Mpc, respectively) are similar to the correlation lengths of 
the reddest and most massive galaxies at similar redshifts (z ~ 2 
and z ^ 1 respectively). 

iMeneux et aD 11006^ find ro values of 2.5 h ^ Mpc for blue 
galaxies and 4.8/i~^ Mpc for red galaxi es at z ~ 0.8 in the Vimos- 
VLT Deep Survey (VVDS). ICoil et alj ( |2007|) find slightly higher 
values in the DEE P2 survey, 3.9 and 5.2h~ ^ Mpc for blue and red 
galaxies at z ~ 1. lAdelberger et all ( l2005b measure ro = 4.5 ± 
0.6,4.2±0.5/i"^Mpcand7 = 1.6 at z = 1.7 and 2.2 for galaxies 
selected through the Lyman-break technique. 

The auto-correlation for QSOs at similar redshifts is less well 
know n: ro = 6.5 ± 1.6h'^ Mpc and 7 = 1.58 ± 0.2 dShen et all 
l2007h: 5.55±0.29 /i~^ Mpc and 7 = 1.63 ± 0.054 on large scales 
dCroom et ^|200 5'. but see this paper for caveats on the measured 
shape of the quasar auto-correlation function). 

Taking the geometric mean of the QSO and galaxy auto- 
correlation amplitudes assuming rg^ = 4.5 and r^'' — 6h~^ Mpc 
gives rQ^ = 5.2h~^ Mpc, which matches our absorber-QSO cross 
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correlation amplitude|f|. Converting ro into tlie alt ernative mea- 
sure o f comoving clustering length Ag, and following lMo & White! 
( l2002h results in QSO dark matter halo masses of > lO^^ '^Mo at 
z ^ 2 and > IO^^Mq at z ~ 1. C IV and Mg II absorption line 
systems are then typically located in galaxies with dark matter halo 
masses of > 1O""^^'^M0. 

8.5 Caveats and Limitations 

The primary limitation in the present analysis is the unknown den- 
sity distribution of absorber clouds in the immediate vicinity of 
the QSO. Our toy model assumes a simple cut-off radius inside of 
which no absorbers can exist, and beyond which power-law clus- 
tering dominates. In a more realistic model, there might be variable 
cut-off radius (7?cut) dependent on the density of the absorbing 
cloud. Such models are beyond the scope of the present analysis. 

We note that a small offset remains between the line-of-sight 
distributions predicted by our model and the the observed one. The 
magnitude of the observed offset is around 200km/s or 2.85 Mpc, 
and may be accounted for by the fact that galaxies in the vicinity 
of the QSO will have a net motion caused by infall into potential 
well of the QSO host halo. This will alter the observed shape and 
peak position of the absorber distribution. A full treatment of infall 
requires detailed cosmological simulations, and will be addressed 
in future work. 

On the observational side, it is clear that a careful correction 
of the QSO redshifts for the effect of the blue-shifting of emis- 
sion lines is warranted in order to constrain the detailed shape of 
the C IV absorber distribution. Follow-up observations in the near 
infra-red of the [O II] and [O III] emission lines in a relatively large 
sample of high-z QSOs would provide a very useful basis for mea- 
suring these effects. The increase in number of QSOs from SDSS 
DR3 to SDSS DR7 will also further enable us to constrain the pre- 
cise shape of the velocity distribution, and the fraction of dense 
absorbers which survive the ionising radiation of the QSO. 



9 SUMMARY 

We have used a cross-correlation analysis of QSO-absorber pairs to 
measure the strength of narrow absorber clustering around QSOs. 
A simple model to convert the 3-D distribution of QSO-absorber 
separations into a line-of-sight distribution in velocity space is pre- 
sented. Our modelling allows us to reach the following conclusions 
for C IV systems: 

• Galaxies in the vicinity of the QSO may contribute as much 
as ~55% of the central spike of absorbers with Aw < 0.01c. The 
remaining ~45% can be directly linked to the presence of the QSO. 

• A high-velocity tail of narrow absorbers is observed, which 
cannot be explained by galaxy clustering. These intrinsic absorbers 
make up ~45% of narrow absorption line systems with 0.01 < 
/\v/c < 0.04. 

• The high-velocity tail is visible in both RQ and RL QSOs. RL 
QSOs show an enhanced number of absorbers at low velocities. 

For Mg II absorbers we find: 

• The spike of "associated" absorbers close to the QSOs in ve- 
locity space is consistent with being caused by galaxy clustering, 

^ The small difference in 7 makes a negligible difference to this result. 



but our analysis cannot rule out the possibility that a considerable 
fraction are caused by absorbers intrinsic to the QSO and its host. 

• The QSO destroys Mg II clouds well beyond the expected 
scale of its own halo, out to at least 800 kpc (comoving). 

• A high-velocity tail of intrinsic Mg II absorbers is detected, at 
a level of 20±5% that observed in the C IV absorbers. 

In the future, the larger absorber samples provided by later 
releases of the SDSS survey data, improved methods for obtain- 
ing reliable QSO redshifts, and investigations of ionisation and line 
width trends with velocity, will contribute substantially to isolat- 
ing the physical processes responsible for QSO outflows detected 
through narrow absorption lines. 
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